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ABSTRACT: We previously identified quantitative trait loci (QTL) for body weight and average daily gain in a common region 
between ADL0198 (chr 1: 171.7 Mb) and ABR0287 (chr 1: 173.4 Mb) on chicken chromosome 1 in an F 2 resource population produced 
by crossing low- and high-growth lines of the Hinai-dori breed. Motilin receptor (MLNR) is a candidate gene affecting growth traits in 
the region. In this study, we genotyped polymorphisms of the MLNR gene and investigated its association with growth traits in a Hinai- 
dori F 2 intercross population. All the exons of the MLNR gene in the parental population were subjected to PCR amplification, 
nucleotide sequenced and haplotypes identified. To distinguish resultant diplotype individuals in the F 2 population, a mismatch 
amplification mutation assay was performed. Three haplotypes (Haplotypes 1-3) were accordingly identified. Six genotypes produced by 
the combination of three haplotypes (Haplotype 1, 2, and 3) were examined in order to identify associations between MLNR haplotypes 
and growth traits. The data showed that Haplotype 1 was superior to Haplotype 2 and 3 in body weight at 10 and 14 weeks of age, 
average daily gain between 4 and 10 weeks, 10 and 14 weeks, and 0 and 14 weeks of age in female in F 2 females. It was concluded that 
MLNR is a useful marker of growth traits and could be used to develop strategies for improving growth traits in the Hinai-dori breed. 
(Key Words: Chicken, Hinai-dori, Growth, Marker Gene, Motilin Receptor) 



INTRODUCTION 

The Hinai-dori is a slow-growing breed of chicken 
native to Akita Prefecture, in northern Honshu Island, Japan. 
Although Hinai-dori meat has a characteristic taste and the 
breed has been used for a long time, it has decreased in 
numbers in recent times owing to the introduction of exotic 
breeds, and for a while was at risk of extinction. The Hinai- 
dori breed has been conserved by hobbyists who belong to 
the Preservation Society (PS) of the Hinai-dori Breed and is 
now mainly used for exhibition purposes. For effective use 
of the breed, selection experiments have been performed at 
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the Akita Prefectural Livestock Experiment Station (LES) 
since 1973 when fertilized eggs were introduced to LES 
from PS with a view towards increasing growth 
performance (Rikimaru et al., 2013). At present, the sex- 
averaged body weight of LES line at 14 weeks of age is 
approximately 0.8 kg heavier than that of PS line (Rikimaru 
et al., 2013). F[ chickens produced by crossing the 
improved LES Hinai-dori sires with Rhode Island Red 
dams have been commercialized as the Hinai-jidori chicken, 
which is one of the most popular high-quality chickens on 
the Japanese market (Rikimaru et al., 2007). 

The identification and utilization of potential candidate 
genes for quantitative trait loci (QTL) with significant 
effects on economically important traits is becoming 
increasingly important in animal breeding programs. There 
are currently two main strategies for detecting QTLs: 
genome scans based on linkage mapping in a cross 
population and association tests using candidate genes 
(Andersson, 2001). In the past decade, QTL mapping for 
chicken growth traits such as body weight has been widely 
studied (Abasht et al., 2006; Hu et al., 2007). Moreover, 
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there have been numerous association studies investigating 
the relationship between growth traits and candidate genes, 
including pre-B-cell colony enhancing factor I (PBEF1) 
(Han et al., 2011), high mobility group AT-hook 2 (HMGA2) 
(Song et al., 2011), insulin-like growth factor I (IGF1) 
(Zhou et al., 2005), pituitary-specific positive transcription 
factor 1 (PIT1) (Nie et al., 2008), and thyroid hormone 
responsive spot 14a (THRSP) (Cao et al., 2007) on 
chromosome 1; insulin-like growth factor binding protein 
(IGFBP) 1 and 3 (Qu et al., 2009) on chromosome 2; 
ornithine decarboxylase (ODC) (Uemoto et al., 2011) on 
chromosome 3; cholecystokinin type A receptor (CCKAR) 
(Rikimaru et al., 2012) on chromosome 4; calpain 3 (Zhang 
et al., 2009) on chromosome 5; IGFBP 2 (Lei et al., 2005; 
Li et al., 2006; Leng et al., 2009) on chromosome 7; growth 
hormone secretogogue receptor (GHSR) (Fang et al., 2010) 
on chromosome 9; and insulin-like growth factor 1 receptor 
(IGF1R) (Lei et al., 2008) on chromosome 10. 

The LES line was developed from PS chickens: hence, 
these individuals possessed genes that influence growth 
traits. To identify QTLs associated with growth traits in the 
Hinai-dori breed, an F 2 resource population produced by 
crossing low- and high-growth lines of the breed was 
analyzed. A significant QTL for body weight at 10 and 14 
weeks of age and average daily gain between 4 and 10 
weeks and between 10 and 14 weeks of age were 
accordingly mapped in a common region between ADL0198 
(chr 1: 171.7 Mb) and ABR0287 (chr 1: 173.4 Mb) on 
chromosome 1 (Rikimaru et al., 2011). Thus, the above- 
mentioned genes, e.g., IGFBP 1-3, ODC, CCKAR, GHSR, 
and IGF1R, are not objects of our interest, since they are 
not located on chromosome 1. Moreover, PBEF1 (chr 1: 
15.2 Mb), HMGA2 (chr 1: 36.1 Mb), IGF1 (chr 1: 57.3 Mb), 
PIT1 (chr 1: 96.2 Mb), THRSP (chr 1: 197.8 Mb) are not 
located between ADL0198 and ABR0287. 

Of all the candidate genes for growth traits on 
chromosome 1, we herein focused our investigation on the 
motilin receptor (MLNR, chr 1: 173.4 Mb) because it was 
located just under the QTL peak adjacent to ABR0287 on 
chromosome 1, and motilin is a key regulatory hormones 
for appetite (Garthwaite, 1985; Rosenfeld and Garthwaite, 
1987), and thus MLNR is a positional and functional 
candidate gene. Our main objective in this study was to 
analyze polymorphism of the MLNR gene and test its 
association with growth traits in Hinai-dori F 2 crossbreds. 

MATERIALS AND METHODS 

Resource population 

PS and LES chickens of the Hinai-dori are shown in 
Figure 1. Fi chickens were produced by crossing three 
cocks from a low-growth line (obtained from PS) with nine 
hens from a high-growth line (obtained from LES). In this 




Figure 1. Hinai-dori. (A) Hinai-dori chickens of the Preservation 
Society (PS) line. Left, male; right, female. (B) Hinai-dori 
chickens of the Akita Prefectural Livestock Experiment Station 
(LES) line. Left, female; right, male. 



cross, one to three hens were randomly selected to mate 
with each male. The F 2 individuals were produced by full- 
sib mating of 17 F[ cocks and 60 F, hens. We accordingly 
obtained 418 F 2 individuals comprising 206 cocks and 212 
hens. The F 2 chickens were hatched on the same day, raised 
in the same chicken house, and fed the same diet ad libitum 
for the duration of the experiment. All animals received 
humane care as outlined in the Guidelines for Proper 
Conduct of Animal Experiments, Science Council of Japan 
(Science Council of Japan, 2006). 

Traits 

Body weight was measured at day 0 (BW 0 d), 4 weeks 
(BW 4 wk), 10 weeks (BW 10 wk), and 14 weeks (BW 14 
wk) of age. Average daily gain between 0 and 4 weeks of 
age (ADG 0-4 wk), between 4 and 10 weeks of age (ADG 
4-10 wk), between 10 and 14 weeks of age (ADG 10-14 
wk), and 0 and 14 weeks of age (ADG 0-14 wk) was 
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Table 1. Polymerase chain reaction primers for the amplification of the 2 exons in the gene encoding the motilin receptor 



Exon 


Set 


Forward primer (5' — » 3') 


Reverse primer (5'— »3') 


Product 
(bp) 


Denaturizing 
temp. (°C) 


Annealing 
temp. (°C) 


1 


1 


TCTTTTTCGTCCTTCGTTCCTT 


CTGAAGTCGGTGCGGTTGT 


704 


99.5 


64.0 


1 


2 


GTCCGACCTGCTCATCCTGCT 


CTCATCCGTCTCACCTCCCCTT 


736 


98.0 


60.0 


2 


1 


GGGGAGACACCCAAGGTCAGGC 


CCAAGCATTTCAGCACACAAAGAGT 


644 


94.0 


62.0 


2 


2 


CCAGGACCTACACAGAGACCAG 


TACAGGGTTGGAGTTAGATGGC 


898 


94.0 


62.0 



calculated from BW at each week of age. 

Identification of MLNR haplotypes 

A draft sequence of the chicken genome (May 2006 
assembly), available on the University of California, Santa 
Cruz (UCSC) Genome Browser (2004) and the Ensembl 
Genome Browser (2004), was used in this study. The 
nucleotide sequences of the two exons of MLNR in the 
parent individuals were determined by polymerase chain 
reaction (PCR) amplification followed by direct sequencing 
to determine nucleotide variance of the gene in the resource 
family. Each genomic DNA of the resource family was 
purified using a SepaGene kit (EIDIA, Tokyo, Japan). Four 
PCR primers were designed to separately amplify the two 
exons of MLNR and the PCRs were performed (Table 1). 
For the PCRs, we used 15-|iL reaction volumes containing 
the following: 10 pmol of each primer for each marker, 200 
^iM of each deoxyribonucleotide triphosphate (dNTP), 1.2 
mM magnesium sulfate (MgSO<0, 0-5 units of KOD-plus 
polymerase (Toyobo, Tokyo, Japan), lx reaction buffer 
provided by the manufacturer, and 10 ng genomic DNA. 
Reactions were performed in a 96 -well plate in an iCycler 
Thermal Cycler (Bio-Rad Laboratories, Hercules, CA, 
USA) using the following conditions: initial denaturation at 
94°C for 2 min; 30 cycles at x°C (x was the denaturizing 
temperature shown in Table 1) for 15 s, at y°C (y was the 
annealing temperature shown in Table 1) for 30 s, at 68 °C 
for 30 s; and a final elongation at 68°C for 4 min 30 s. The 
PCR products were purified using High Pure 96 UF 
Cleanup Plates (Roche Diagnostics GmbH, Mannheim, 
Germany) and subjected to cycle sequencing using a 
BigDye Terminator Cycle Sequencing FS Ready Reaction 
kit (Perkin-Elmer Applied Biosystems, Foster City, CA, 
USA) and the same primers as used for the PCR 
amplification. Sequencing was performed in both directions 
and each sequencing product was run in an ABI 3130 
Genetic Analyzer (Perkin-Elmer Applied Biosystems). The 
MLNR gene haplotypes in the F 2 intercross population were 
then identified. 

Statistical analyses 

A mixed-inheritance animal model was used to evaluate 
the effects of haplotypes on traits. The ld_fix option of the 
Qxpak software (Perez-Enciso and Misztal, 2004) was used 
for the MLNR haplotypes. For statistical models, the 



additive effect of each allele and the sex effect were used as 
fixed effects. In this analysis, the infinitesimal genetic effect 
was included and treated as a random effect with covariance 
matrix Acr^ (where A is the numerator relationship 

matrix). Likelihood ratio tests were performed by not 
including the additive effects in the model, and nominal p- 
values were obtained by assuming a x 2 distribution of the 
likelihood ratio test. 

Percentage of haplotype variance explained by the 
model was calculated as 

Variance percentage 

= 100x(R_variance-F_variance)/R_variance 

where R_variance is the residual and polygenic 
variances from the reduced model, omitting haplotype 
effect but including sex and polygenic effects, and 
F_variance is the residual variance from the full model, 
including haplotype, sex, and polygenic effects. 

RESULTS 

On sequencing of the MLNR gene, we found GGGGC 
insertion polymorphisms located 54 bp upstream of the 5'- 
untranslated region and the sequenced region differed by 6 
single nucleotide polymorphisms (SNPs) from the 
corresponding region in the draft sequence of the chicken 
genome. Resultantly, three MLNR haplotypes (Haplotypes 
1-3) were identified in the parent individuals (Table 2). The 
nucleotide sequences of the three haplotypes have been 



Table 2. Genotypes of parent individuals of the F 2 populations 


Individuals 


No. 


Genotype 


Parent male 


1 


2/2 




2 


2/2 




3 


2/3 


Parent female 


1 


1/2 




2 


1/2 




3 


2/2 




4 


2/2 




5 


2/2 




6 


1/2 




7 


1/2 




8 


2/2 




9 


1/1 



Takahashi et al. (2014) Asian Australas. J. Anim. Sci. 27:316-323 31 9 



Table 3 Polymerase chain reaction (PCR) primers for the mismatch amplification mutation assay (MAMA) 



Primer sets 


Primers (5'->3') 


PCR product 
(bp) 


Target position 
(AB703068) 


1 


PCR amplification 
Haplotype 

2 


3 


MLNR-1008-13 


CAGCGTCCTCTACGGCTTCA 


217 


gl008. C>T 


+ 




+ 




GCCGGCGGTCCGTGGtCa 












MLNR-1008-2 


CAGCGTCCTCTACGGCTTCA 


217 


gl008. C>T 




+ 






GCCGGCGGTCCGTGGtCg 












MLNR- 1029-1 


GCGTCCTCTACGGCTTCATC 


238 


gl029. A>G 


+ 








AGTCTGTGGAATGGCGCgTc 












MLNR- 1029-23 


GCGTCCTCTACGGCTTCATC 


238 


gl029. A>G 




+ 


+ 




AGTCTGTGGAATGGCGCgTt 













Bases shown in lower case with an underline represent induced mismatches. Bases shown in lower case at the 3 '-end represent target single nucleotide 
polymorphisms (SNPs). The haplotypes that can or cannot be amplified by PCR for each primer set are shown as '+' or respectively. For the PCRs, 
we used 15-pL reaction volumes containing the following: 3 pmol of each primer for each marker, 200 pM of each dNTP, 0.5 units of Paq5000DNA 
Polymerase (Agilent Technologies, La Jolla, CA, USA), lx reaction buffer (containing 2 mM MgCl 2 ) provided by the manufacturer, and 10 ng genomic 
DNA, Reactions were performed in a 96-well plate in an iCycler Thermal Cycler (Bio-Rad Laboratories, Hercules, CA, USA) using the following 
conditions: initial denaturation at 95°C for 2 min; 30 cycles at 94°C for 30 s, at 62°C for 30 s, at 72°C for 1 min (MLNR-1008-13 and -2), or 35 cycles at 
94°C for 30 s, at 64°C for 30 s, at 72°C for 1 min (MLNR- 1029-1 and -23); and a final elongation at 72°C for 7 min. The PCR products were 
electrophoresed on a 2.0% agarose gel with lx Tris-acetate EDTA (TAE) buffer and stained with ethidium bromide. The combination of these results 
enabled us to identify diplotypes in the F 2 population. 



registered in the DNA data bank of Japan (DDBJ) with the 
following accession numbers: AB703068, AB703069, and 
AB703070. The results indicated that these haplotypes 
might be segregated and individuals with six different 
genotypes (1/1, /2/2, 3/3, 1/2, 1/3, and 2/3) might appear in 
the F 2 population. Therefore, a mismatch amplification 
mutation assay (MAMA) PCR protocol was developed that 
detects the 6 genotypes described by Cha et al. (1992). We 
designed four PCR primers to distinguish MLNR haplotypes 
and PCR and genotyping were performed as described in 
Table 3 and Figure 2. 

We observed an association between three MLNR 
haplotypes (Haplotype 1, 2, and 3) and growth traits (Table 
4). There was no association between these haplotypes and 
growth traits in the overall F 2 population. When data were 

1008-13 1008-2 

I ) I \ 

M 1 2 3 4 



500 bp^ 
400 bp- 

200 bp 1 

100 bp- 




analyzed separately by sex, a significant (p<0.05) 
association was found between these haplotypes and BW 10 
wk, BW 14 wk, ADG 4-10 wk, ADG 10-14 wk, and ADG 
0-14 wk in F 2 females, whereas there was no association 
between these haplotypes and growth traits in F 2 males. 

DISCUSSION 

Motilin, a 22 amino acid polypeptide first discovered in 
the mucosa of upper small intestine of the hog (Brown et al., 
1971), is predominantly recognized for its prokinetic effects 
on gastrointestinal motility and emptying (Poitras et al., 
1994). Motilin has been identified in many mammals (Bond 
et al., 1988; Banfield et al., 1992; Xu et al., 2001; 
Strausberg et al., 2002; Xu e al., 2003; Ohshiro et al., 2008) 

1029-1 1029-23 



500 bp 
400 bp 

200 bp 
100 bp 




Figure 2. Typical patterns of 2% agarose gel electrophoresis of MAMA results. (A) MAMA results at the SNP site of AB703068: g.1008 
C>T. MLNR-1008-13 and MLNR -1008-2 primer sets were used in lane 1 and 2, and in lane 3 and 4, respectively. Lane M, DNA size 
marker; lane 1 and 3, patterns of Haplotype 2; lane 2 and 4, patterns of Haplotype 1 and 3. (B) MAMA results at the SNP site of 
AB703068: g.1029 A>G. MLNR-1029-1 and MLNR-1029-23 primer sets were used in lane 5 and 6, and in lane 7 and 8, respectively. 
Lane M, DNA size marker; lane 5 and 7; pattern of Haplotype 2 and 3; lane 6 and 8; patterns of Haplotype 1. 
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Table 4. Phenotypic values of growth traits and effects of Motilin receptor haplotypes on the traits in the F 2 population 



Trait 


n 


Phenotypic values 


LRT 


p value 


Haplotype 1 


Haplotype 2 


Haplotype 3 


Variance 


Mean±SD 


meaniSE 


meaniSE 


meaniSE 


(%) 


Total 


















BW 0 d (g) 


418 


34.4±2.4 


2.7 


0.255 


_ 


_ 


_ 


_ 


BW 4 wk (g) 


417 1 


231.1+38.3 


0.2 


0.888 










BW 10 wk (g) 


418 


960.61163.5 


4.3 


0.116 










BW 14 wk (g) 


418 


1,466.9+260.7 


5.2 


0.075 


_ 


_ 


_ 


_ 


ADG 0-4 wk (g/d) 


417 1 


5.2+1.0 


0.3 


0.856 










ADG 4- 10 wk (g/d) 


417 1 


17.4+3.3 


5.1 


0.077 










ADG 10-14 wk (g/d) 


418 


18.1+4.4 


3.6 


0.166 










ADG 0-14 wk (g/d) 


418 


14.6+2.7 


5.2 


0.073 










Male 


















BW 0 d (g) 


206 


34.6±2.5 


5.0 


0.041 


_ 


_ 


_ 


_ 


BW 4 wk (g) 


205 1 


243.6±38.3 


0.3 


0.797 










BW 10 wk (g) 


206 


1,053.4+144.8 


2.2 


0.536 










BW 14 wk (g) 


206 


1,658.8+199.5 


1.5 


0.670 










ADG 0-4 wk (g/d) 


205 1 


5.5+1.0 


0.3 


0.835 


_ 


_ 


_ 


_ 


ADG 4-10 wk (g/d) 


205 1 


19.3+2.9 


2.8 


0.485 


_ 


_ 


_ 


_ 


ADG 10-14wk(g/d) 


206 


21.6+2.9 


0.1 


0.970 










ADG 0-14 wk (g/d) 


206 


16.6+2.0 


0.8 


0.672 


_ 


_ 


_ 


_ 


Female 


















BW 0 d (g) 


211 


34.2±2.3 


0.7 


0.915 


- 


- 


- 


- 


BW 4 wk (g) 


211 


219.0±34.3 


3.2 


0.374 










BW 10 wk (g) 


211 


870.81125.9 


8.6 


0.021* 


462.9120.2 


426.1115.5 


394.4123.8 


7.5 


BW 14 wk (g) 


211 


1,280.6+157.4 


12.0 


0.002* 


692.6125.6 


629.8119.5 


582.5130.3 


11.0 


ADG 0-4 wk (g/d) 


211 


4.910.9 


3.0 


0.375 










ADG 4-10wk(g/d) 


211 


15.512.5 


8.9 


0.041* 


8.210.4 


7.610.3 


7.010.5 


5.7 


ADG 10-14 wk (g/d) 


211 


14.612.3 


7.5 


0.007* 


8.210.4 


7.310.3 


6.710.5 


5.5 


ADG 0-14 wk (g/d) 


211 


12.711.6 


12.3 


0.002* 


6.910.3 


6.310.2 


5.810.3 


11.4 



BW 0 d, BW 4 wk, BW 10 wk, and BW 14 wk: body weight measured at day 0, 4 weeks, 10 weeks, and 14 weeks of age. 

ADG 0-4 wk, ADG 4-10 wk, ADG 10-14 wk, and ADG 0-14 wk: average daily gain between 0 and 4 weeks, between 4 and 10 weeks, between 10 
and 14 weeks, and between 0 and 14 weeks. 

SD = Standard deviation. SE = Standard error. LRT = Log likelihood ratio test statistics. 

1 Since BW 4 wk has a missing datum, ADG 0-4 wk and 4-10 wk have missing data. * p<0.05. 



and chickens (DeClercq et al., 1996). Motilin binding sites 
are present in the gastrointestinal tract and central nervous 
system, along with a G-protein-coupled receptor, GPR38, 
which has been identified as MLNR (Feighner et al., 1999). 
Initially, GPR38/MLNR was identified as an orphan 
receptor belonging to the GHSR family. This family 
contains receptors for ghrelin (Kojima et al., 1999), 
neuromedin U (Howard et al., 2000), and neurotensin 
(Vincent et al., 1999). Of the receptors, the GHSR is most 
homologous to the MLNR; the chicken forms share 53% 
identical amino acids (Yamamoto et al., 2008). 

It has been shown that motilin and ghrelin have 
contractile activities in the gastrointestinal tract and their 
receptors are presented in the gastrointestinal tract of 
chickens (Kitazawa et al., 1997; Tanaka et al., 2003), 
although detail distributions of MLNR and GHSR have not 
been reported. DeClercq et al. (1996) reported that motilin 



was most abundant in the mucosa of the colon, followed by 
the jejunum, duodenum, and ileum; low motilin contents 
were found in the proventriculus and gizzard. Neglia et al. 
(2004) reported that the largest numbers of ghrelin-positive 
cells were located at the base of lobuli of the proventriculus 
gland, along villi of the intestines and in crypts of the 
duodenum, while lower numbers of ghrelin-immunostained 
cells were located in crypts of jejunum and ileum and only 
few ghrelin-immunostained cells were detected at the base 
of crypts of the large intestine. Kitazawa et al. (2007) 
reported that chicken ghrelin caused contraction of chicken 
gastrointestinal smooth muscles in a region dependent 
manner (crop>colon>esophagus>proventriculus>jejunum> 
ileum>duodenum), and chicken motilin caused contraction 
of chicken gastrointestinal tract in a region dependent 
manner (duodenum>j ej unum>ileum>colon>pro ventriculus 
>crop>esophagus). These data suggest that ghrelin mainly 
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regulates contractions of the upper (crop, esophagus, and 
proventriculus) gastrointestinal tract, motilin principally 
regulates small intestinal motility in chicken. Since it has 
been reported that motilin can stimulate the GHSR but 
ghrelin does not activate MLNR in rabbits (Dass et al., 
2003) and chickens (Yamamoto et al., 2008), motilin 
abundantly presented in colon may activate GHSR and lead 
gastrointestinal mortality in colon. Taken together, these 
reports suggest that motilin and MLNR play central roles in 
the regulation of gastrointestinal motility and gastric 
emptying with ghrelin and GHSR in chicken. 

The association between haplotypes and growth traits 
cannot be explained by the protein sequences, since the 
SNPs detected in the coding regions of MLNR do not 
induce missense mutations, and thus MLNR proteins are 
functionally normal. Therefore, at present, it is unknown 
whether MLNR haplotypes have any functional significance. 
The associations detected in this study may be produced by 
linkage disequilibrium between these haplotypes and 
another linked gene directly involved in the regulation of 
these traits. Recently, Zhang et al. (2011) proposed 
positional candidate genes for growth traits in the same 
region on chromosome 1, namely, retinoblastoma 1 (RBI) 
(chr 1: 173.0 Mb); lysophosphatidic acid receptor 6 (chr 1: 
173.0 Mb); fibronectin type III domain containing 3A (chr 
1: 173.3 Mb); MLNR (chr 1: 173.4 Mb); and calcium 
binding protein 39-like (chr 1: 173.5 Mb). Of the genes, 
Zhang et al. (2011) suggested that two SNPs of RBI gene 
might be quantitative trait nucleotides for BW; however 
they have not proposed any hypothetical mechanisms why 
the SNPs affect BW. 

One likely possibility to explain the association between 
MLNR haplotypes and growth traits is that presence of 
GGGGC insertion polymorphism in the 5 '-upstream region 
of MLNR, corresponding to the predicted Spl binding site, 
may affect MLNR gene expression. Spl is a transcription 
factor that plays an important role in many cellular 
activities such as metabolism, cell growth, differentiation, 
angiogenesis and apoptosis (Deniaud et al., 2009). Spl 
contains three Cys2His2-type zinc finger motifs and it is 
generally accepted that its mode of action is mediated 
mainly by binding to the decanucleotide consensus 
sequence 5'-(G/T)GGGCGG(G/A)(G/A)(C/T)-3' in double- 
stranded DNA (dsDNA) (Song et al., 2001). In the 
predicted Spl binding site of MLNR gene, Haplotype 1, 2, 
and 3 have (GGGGC) 3 , (GGGGC) 4 , and (GGGGC) 5 , 
respectively, while the draft sequence of the chicken 
genome has (GGGGC) 2 . Since it was reported that Spl- 
tandem repeat polymorphisms affect binding capabilities of 
Spl and transcriptional efficiencies (In et al., 1997; 
Sigurdsson et al., 2008), the repeat polymorphisms in 
MLNR gene promoter may affect MNLR gene expression. 
Our data suggest that shorter GGGGC-repeats shows better 



growth performance rather than longer ones. Therefore, 
further studies are needed to elucidate the differences in 
transcriptional and translational efficiencies caused by 
haplotypes. 

It is also unclear the reason that MLNR haplotypes are 
associated with gender differences in growth traits. Motilin 
and MLNR genes in mice and rats exist as pseudogenes and 
are thought to genetically knockdown; however, it is still 
debated whether or not native MLNR is expressed (He et al., 
2010). Recently, Miegueu et al. (2011) demonstrated that 
motilin stimulates preadipocyte proliferation and 
differentiation and adipocyte lipid storage in murine 3T3- 
Ll and primary rat adipocytes. Previous reports indicated 
that acylated and unacylated ghrelin might influence 
preadipocyte proliferation (Zhang et al., 2004) and could 
enhance preadipocyte differentiation (Miegueu et al., 2011), 
thus, motilin may act via GHSR in the murine and rat 
adipocytes. Zigman (Zigman et al., 2005) reported that 
female GHSR-null mice were significantly lighter than 
wild -type controls between 12 and 19 weeks of age, while 
male GHSR-null and wild-type mice did not differ in body 
weight at any weeks of age. After 19 weeks of age, female 
GHSR-null mice had 35.6% less body fat than wild-type 
mice. Taken together, these data may imply that MLNR 
plays a key role in fat deposition and body weight 
especially in females; however, further studies needed to 
clarify the reason how MLNR affect fat traits in chicken. 

In conclusion, we observed significant associations 
between MLNR haplotypes and growth traits in an F 2 
population produced by crossing low- and high-growth 
lines of the Hinai-dori breed of chicken. Our data suggest 
that haplotype-assisted selection is a potentially useful 
approach in the Hinai-dori breed. In the future, the 
associations between haplotypes and growth traits should be 
studied in other chicken lines and breeds. 
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